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ABSTRACT 


The  detonation  velocity  for  a  granular  mixture  of  potassium  per¬ 
chlorate  and  aluminum  confined  in  a  Lucite  tube  and  an  ideal  detonation 
velocity  calculated  by  the  RUBY  Computer  Program  have  been  used  to¬ 
gether  with  the  observed  reaction  zone  profile  to  estimate  the  position 
of  the  Chapman-Jouguet  plane  with  respect  to  the  plane  at  which  reaction 
commences.  This  was  done  by  extending  Jones  equation  for  effect  of 
radial  expansion  on  the  detonation  velocity  to  include  the  effect  of 
covolume  a  function  of  volume.  It  was  found  that  the  CJ  plane  is  about 
0,9  cm.  behind  the  plane  at  which  reaction  commences. 
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INTRODUCTION 


A  cylindrical  charge  of  explosive  detonates  at  a  lower  velocity 
than  the  ideal  detonation  velocity  since  the  latter  corresponds  to  a 
one  dimensional  process.  The  ;mlocity  decrease  is  governed  by  -  the 
cross-sectional  area  increase  at  the  Chapman- Jouguet  (CJ)  plane  which 
depends  on  the  confinement  of  the  charge,  the  equation  of  state  of  the 
products,  and  the  fluid  dynamical  equations,  Jones  (Ref,  l)  has  treated 
the  case  where  a  constant  covolume  equation  of  state  may  be  assumed. 

The  granular  explosive,  considered  in  this  report,  has  a  long  re¬ 
action  zone,  and  doubles  cross  section  in  about  a  centimeter  of  length. 
Since  a  constant  covolume  can  not  be  assumed  for  this  case,  the  approach 
of  Jones  has  been  extended  for  covolume  a  function  of  volume.  The 
equations  developed  permit  determining  area  increase  at  the  CJ  plane  as 
a  function  of  the  ratio  of  observed  to  ideal  detonation  velocity.  This 
was  dene,  using  for  the  latter  a  theoretical  value  calculated  with  the 
RUBY  Computer  Program  (Ref,  2),  The  area  increase  so  derived  was  com¬ 
pared  with  the  reaction  zone  profile  to  locate  the  CJ  plane. 

The  granular  explosive  under  study  (Ref,  3-5)  is  a  mixture  of  po¬ 
tassium  perchlorate  (22  micron)  and  atomized  aluminum  (15  micron)  in 
the  weight  ratio  60/40,  For  convenience,  this  mixture  is  hereafter 
referred  to  as  60/40,  The  powder  is  loose-loaded  (loading  density 
1*5  g/cm3)  into  a  3/8  inch-diameter  axial  hole  in  a  20  inch  long,  3  inch 
diameter  Incite  tube.  Initiation  is  by  a  confined  0,36  gram  tetryl 
pellet,  A  stable  propagation  profile  at  a  constant  velocity  of  about 
900  m/aec  is  observer  for  the  latter  part  of  the  tube,  A  typical  stable 
profile  is  shown  in  Figure  1,  Photographs  showing  the  development  and 
progression  of  the  detonation  are  in  References  3  and  5. 
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EQUATIONS 


Jones  (Ref,  l)  derives  the  following  system  of  equations  in  terms 
of  the  non-dimensional  coordinates  y,  z,  and  other  symbols  defined  be¬ 
low.  Details  are  available  in  the  cited  reference.  The  introduction 
of  the  additional  variable  r  means  the  system  is  short  one  equation 
for  a  unique  solution,  Jones  formulation  recasts  the  basic  equations 
so  that  one  variable  y  can  be  considered  a  constant,  thus  permitting 
a  particular  solution. 


Continuity 


Chapman-Jouguet  Condition 


(4) 
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where 


r  denotes  the  ratio  of  the  radius  of  a  stream  tube  at  any  point 
to  that  at  the  front,  r2  is  a  measure  of  radial  area  exp«4nsion. 

p  denotes  the  pressure  in  excess  of  one  atmosphere. 

U  is  the  velocity  of  the  detonation  wave, 
is  the  loading  density  of  the  explosive. 

P  is  the  density  of  the  fluid  at  any  point  beyond  the  shock  wave 

^  front. 

q  is  the  particle  (fluid)  velocity  behind  the  front,  relative  to 
the  detonation  wave  at  rest, 

f  is  a  dimensionless  function  of  the  non-dimensional  coordinates 
y,  z  and  the  degree  of  reaction  €  . 

E  is  the  internal  energy  per  unit  mass  relative  to  the  initial 
site. 


()  is  the  chemical  energy  liberated  per  unit  :/>ass  by  the  reaction, 
k  is  a  parameter  in  the  equation  of  state  related  to  specific  heat, 

b  is  a  parameter  in  the  equation  of  state  related  to  covolume. 


The  subscript  zero  refers  to  conditions  at  the  beginning  of  the  re¬ 
action  zone,  that  is  behind  the  shock.  The  subscript  one  refers  to 


conditions  at  the  Chapssan— Jouguet  plane.  No  subscript  indicates  any 
point  between  these.  2  =:  /b / 

y  is  a  mean  value  of  y  through  the  reaction  zone  defined  in  Eq,  11 
which  follows  later.  It  is  a  new  variable  which  must  be  approxi¬ 


mated  to  obtain  a  solution. 


Uq  is  the  detonation  velocity  when  there  is  no  expansion  of  the 
reaction  zone. 


Jones  treats  k  and  b  as  constants  and  shows  that  the  changes  in 
detonation  velocity  ir  terms  of  radial  expansion  at  the  CJ  plane  are 


This  formula  could  not  be  used  because  the  results  calculated  with  it 
contradicted  the  assumption  that  the  covolume  was  constant.  The  equation 
of  state  was  then  modified  to  treat  b  (related  to  covolume)  as  a  linear 
function  of  y  (the  normalized  density).  This  leads  to  a  modification  of 
Eq.  6  (to  Eq.  9)  as  follows:  Inserting  Eq.  5  in  4  and  then  inserting 
Eq.  2a  in  the  result  leads  to  Eq.  7  below.  Eq.  8  follows  by  using  Eq.  7 
in  Eq.  2a.  b*  is  the  derivative  of  b  with  respect  to  y.  Eqs,  7  and  8 
are  the  analogues  of  Jones  Eqs.  17  and  18. 

(i+  =  A(i  -  i^O'^V)  (8) 


If  we  substitute  Eqs.  5i  7»  and  8  into  Eq.  3a  we  shall  obtain  a 
relation  expressing  U  in  terms  of  k,  b,  b*  ,  r  and  y.  Ely  setting 
U  =  U  for  r  =  1  and  eliminating  between  the  original  and  modified 
equation  we  arrive  at  the  desired  result  for  variable  covolume,  Eq.  9 
below. 


/»-  .5 

ir) 


where 
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-  Z  -f  6  p 
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c 
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-  ^4^  +  4  ^  ~  J 
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(t  +  Zk-b'M) 


(1+34  - 
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Equation  9  reduces  to  Eq.  6  (Jones  equation)  for  ij  =  0  and  b  = 
(covolurae  for  r  =  1,  considered  constant).  Note  also  that  for  r^  =  1, 
^  a  1  that  one  obtains  the  expected  U  =  U^. 

RESULTS 

Equation  9  shows  how  the  detonation  velocity  (U)  depends  on  the 
radial  expansion  at  the  CJ  plane  (r^).  Conversely,  using  the  observed 
U,  and  providing  the  other  quantities  in  Eq.  9i  one  can  find  r^  and 
locate  the  CJ  plane  from  the  observed  reaction  zone  profile.  The  ac¬ 
curacy  of  the  result  of  this  calculation  is  limited  by  the  variability 
in  the  experimental  results  and  the  assumptions  necessary  to  specify 
the  other  quantities  in  the  equation.  In  the  paragraphs  that  follow, 
numerical  values  to  be  used  in  Eq.  9  are  developed. 

Since  b  =  a  ,  where  a  is  the  covolume  per  unit  mass,  and  A  is 

the  loading  density  (1.5  g/cm^),  it  follows  that  b  is  the  covoluroe  in 
units  ot  1/ A  .  To  calculate  values  of  b,  the  RUBY  Computer  Program 
was  used  to  calculate  the  CJ  properties  and  those  along  an  isentrope 
for  the  mixture  under  consideration.  Since  the  program  uses  a 
Kistiakowsky-Wilson  equation  of  state 

JL  =  /  + 

n??r 

and  the  co volume  form  is  P(V  -  a)  =  nRT, 

the  gaseous  covolume  is  represented  in  RUBY  notation  by 

^  "  P  / 

(7 

To  this  was  added  the  volume  per  gram  of  explosive  occupied  by  the  two 
solid  products,  A1  and  AigO^.  The  covoluraes  obtained  are  tabulated 
below  versus  pressure  and  the  RUBY  density  ratij:  RHOLO  which  is  the 
same  as  the  y  in  Jones  notation. 
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RHOLO(y) 

b 

0.7080 

0.6419 

lea 

1.1757 

0.9933 

Ico 

2.0554 

1.5433 

These  values  were  found  to  be  fitted  very  well  by  the  straight  line 
equation 

Thus,  the  value  (1.5)"  =  0.666?  is  appropriate  for  the  derivative  b  . 

The  quantity,  k  =  in  the  equation  of  state  is  the  ratio  of  the 
specific  heat  from  E  =  C^T  to  the  gas  constant  from  p(V  -  a)  =  RT.  The 
value  of  k  was  calculated  for  those  explosive  products  predicted  by  the 
RUBY  program  at  the  CJ  point.  Solid  products  were  included  in  the 
weighted  average  for  the  specific  heat  calculation,  but  not  in  that  for 
the  gas  constant.  The  value  k  =  7.935  resulted.  No  significant  change 
from  this  value  occurred  for  the  products  corresponding  to  a  pressure 
drop  along  the  isentrope  to  l/lO  the  CJ  value.  Thus,  k  was  treated  as 
a  constant  in  the  derivation  of  Bq.  9. 


The  quantity  y  is  defined  in  the  Jones  derivation  as 


It  is  an  average  of  the  non-dimensionalized  specific  volume  over  an  in¬ 
verse,  cross-sectional  area  squared,  coordinate.  The  integration  is 
from  the  point  behind  the  initial  shock  of  the  detonation  zone  where 
radial  expansion  (r  «  l)  begins  to  the  CJ  plane  (r  =  r^).  Vftien  there 
is  no  radial  expansion  then  the  density  at  the  CJ  plane  is  found  from 
the  RUBY  calculation  to  be  2.115  g/cra^.  The  density  behind  the  shock  ( ) 
follows  from  combining  an  estimated  shock  Hugoniot  of  the  powder  with 
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'■i 

the  observed  detonation  velocity.  It  is  about  ?.4  g/cra''.  If  one 

assumes  that,  at  the  CJ  plane,  when  lateral  expansion  has  occurred  the 

2 

value  of  =  2,115  decreased  inversely  as  r^  and  then  uses  a  linear 
average  of  densities,  then  y  is  approximated  by 


An  alternate  approach  is  to  approximate  y  in  Eq,  11  by  having  the  density 

2 

decrease  from  the  value  at  r  =  1  inversely  proportional  to  r  ,  Inte¬ 
gration  of  Eq.  11  then  leads  to 


^  ft*  1 

r  < 

where  ^  5  Y 

Numerical  values  for  these  two  assumptions  are 

listed  below. 

r^  (area  ratio/ 

1.1 

1.5 

2.0 

y  from  Eq,  12 

0.6942 

0.7876 

0.8679 

y  from  Eq,  I3 

0.6548 

0.7500 

0.8333 

(13) 


In  summary,  Eq.  9#  using  the  numerical  values  and  relations  for 
y  and  b  from  the  previous  paragraphs,  now  gives  in  terms  of  r^. 

Further,  from  Eq,  8,  we  can  find  the  corresponding  pressure  at  the  CJ 
plane  and  form  the  ratio  to  that  for  r^  =  1.  The  values  obtained  are 
shown  below  together  with  the  value  of  U  that  corresponds  to  a  theoretical 
ideal  detonation  velocity  (from  RUBY)  =  4966  m/sec. 
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1.5 

2.0 

2.5 

u  /u 

o' 

1.93 

2.69 

3.38 

U(m/sec) 

2569 

1848 

1468 

Pn/Pi 

0.294 

■  0.155 

0.098 

3.0 

3.5 

4.0 

4.5 

4.07 

4.71 

5.35 

5.99 

1221 

1055 

928 

829 

0.067 

0.050 

0.039 

0.031 

It  follows  that  for  an  experimentally  observed  detonation  velocity 
of  928  m/sec  the  CJ  plane  is  located  where  r^  is  4,  or  r^  =  2  and  the 
pressure  drops  from  the  RUBY  value  for  r^  =  l(pj^  =  107  kilobars)  to 
0,039  of  it,  that  is,  to  4.2  kilobars. 


Designate  by  the  distance  measured  away  from  the  front,  starting 
from  the  point  where  the  area  expansion  of  the  reaction  zone  starts. 

Call  the  observed  angle  of  this  expansion,  Q  and  the  cylindrical  radial 
distance  from  the  axis,  d.  Then  for  no  expansion,  d  0.5  cm  which  is 
half  the  diameter  of  the  explosive  column.  For  the  CJ  plane  it  follows 


that 


d--  0,5 


oL 


{o.sY' 


whence 


0-5 


(14) 


To  use  Bq,  14,  the  value  of  O  ,  must  be  obtained  from  the  observed 
reaction  zone  profiles.  Figure  2  shows  the  variation  in  profiles  ob¬ 
served  for  different  loading  conditions.  For  Q  -  JO  degrees  Eq,  14 
gives  =  0.87  cm  which  at  928  m/sec  corresponds  to  9.3  microsec. 
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DISCUSSION 


The  objective  of  the  calculations  presented  herein  was  only  to 
obtain  an  estimate  of  the  location  of  the  Chapman-Jouguet  plane  using 
the  fact  that  the  observed  detonation  velocity  was  one-fifth  the  theo¬ 
retical  ideal  detonation  velocity.  It  was  recognized  that  the  assumptions 
(discussed  later)  necessary  to  obtain  a  solution  only  supported  the 
conclusion  that  the  distance  of  the  CJ  plane,  from  the  point  where  reaction 
commences  is  in  the  range  of  a  fraction  of  a  centimeter  to  a  few  centi¬ 
meters,  (0.8?  cm  was  calculated  corresponding  to  9»3  microsec  for  a 
velocity  of  928  m/sec). 

It  takes  a  minimum  of  15  microsec  to  vaporize  a  15  micron  aluminum 
particle  under  the  detonation  zone  conditions.  (It  is  planned  to  treat 
the  kinetics  of  aluminum  particle  vaporization  in  a  separate  report.) 

Thus,  our  calculation  suggests  incomplete  combustion  of  the  aluminum 
at  the  CJ  plane,  and  supports  the  idea  that  physical  kinetics  is  a  signi¬ 
ficant  factor  in  the  detonation  of  a  granular  explosive. 

Constraints  on  available  fuel  and  associated  chemical  energy  release 
were  not  included  in  the  RUBY  calculation  of  the  ideal  detonation  velocity. 
To  determine  the  effect  of  such  constraints,  the  RUBY  calculation  was 
repeated  with  33  and  675S  of  the  aluminum  available  in  the  explosive  made 
to  appear  as  aluminum  in  the  products.  The  explosive  energy  was  altered 
correspondingly.  Although  numerical  results  were  changed,  the  position 
of  the  CJ  plane  was  still  in  the  range  of  a  fraction  of  a  centimeter  to 
a  few  centimeters. 

The  validity  of  the  result  depends  on  the  assumptions  made  and  the 
distributions  associated,  with  the  numeri.cal  values  used  in  the  calculation. 
The  quantities  involved  are  (U/U^),  ©  »  y.  b,  b^  ,  and  k.  U  and  0  have 
been  observed  in  over  fifty  experiments.  Some  of  these  involved  slightly 
different  mixtures.  It  appears  safe  to  say  that  6  falls  between  15  and 
45  degrees  (30  degrees  was  used)  and  that  U^  is  between  800  and  1100  ra/sec 
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(900  m/sec  was  used).  U  has  been  calculated  for  many  variations  of  input 

conditions,  U  for  conditions  approximating  the  experiments  and  physical 

kinetics  ranges  from  4100  to  5800  s/sec  (5000  m/sec  was  used).  Hence, 

U/U  must  fall  between  4  and  7  (5  was  used)  and  (9  between  15  and  45 
o 

degrees  (30  degrees  was  used). 

The  value  of  k  =  was  a  very  slowly  changing  function  of  the 
RUBY  explosive  products.  The  key  assumption  believed  correct,  is  that 
the  condensed  products  contribute  to  specific  heat,  but  not  to  the  gas 
constant. 

The  value  of  b  was  first  taken  as  constant  and  Jones  formula  used 
directly.  The  result  was  not  much  different  than  that  presented  herein 
for  b  a  function  of  colume.  However,  the  covolume  calculated  for  the 
mixture  that  would  exist  at  the  CJ  plane  was  double  the  constant  value 
used  in  the  calculation.  It  is  for  this  reason  that  the  calculation  was 
redone  for  variable  covolume.  However,  the  fact  that  the  results  are 
similar  for  the  two  calculations  indicates  that  b  and  I5  variability  is 
not  a  controlling  factor. 

The  key  y  assumption  was  that  the  density  that  exists  at  the  CJ 
plane  decreased  inversely  as  the  area  ratio  increased  there.  This  ignores 
(for  y  only)  the  effect  on  density  of  changes  in  particle  velocity 
occasioned  by  the  area  increase.  It  was  found  that  if  y  was  altered  by 
+  50^  in  the  equations  with  the  other  quantities  in  their  normal  range 
that  the  CJ  location  could  be  shifted  in  the  range  from  2  times  to  l/3 
the  computed  value. 

Combining  all  the  considerations  of  the  above  paragraphs  it  was 
decided  that  one  could  only  conclude  that  the  position  of  the  CJ  plane 
relative  to  the  plane  at  which  reaction  commences  was  in  the  range  of  a 
fraction  of  a  centimeter  to  a  few  centimeters.  This  is  adequate  to  show 
that  this  granular  medium  confined  in  a  Incite  tube  has  a  detonation  zone 
two  magnitudes  greater  than  that  for  a  cast  organic  explosive.  The  long 
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detonation  zone  and  information  on  fuel  vaporization  rates  and  observation 
of  afterburning  show  that  physical  kinetics  is  a  strong  factor  in  the 
detonation. 
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Fig.  1  lypical  stable  profile  of  the  detonation  of  a  granular  mixture  of  potassium  perchlorate 
and  aluminum  confined  in  a  Lucite  tube.  The  detonation  velocity  for  this  experiment 
(No.  1024)  was  900  m/sec. 
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Stable  reaction  zone  profiles  observed  in  different  experljnents  on  the  detonation 
of  a  granular  mixture  of  potassium  perchlorate  and  aluminum  confined  in  a  Lucite 
tube.  Numbers  shown  designate  the  experiments,  not  the  detonation  velocities. 


